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A linear stability analysis of the free surface of a horizontally unbounded ferrofluid layer of ar- 
bitrary depth subjected to vertical vibrations and a horizontal magnetic field is performed. A 
nonmonotonic dependence of the stability threshold on the magnetic field is found at high frequen- 
cies of the vibrations. The reasons of the decrease of the critical acceleration amplitude caused by a 
horizontal magnetic field are discussed. It is revealed that the magnetic field can be used to select 
the first unstable pattern of Faraday waves. In particular, a rhombic pattern as a superposition of 
two different oblique rolls can occur. A scaling law is presented which maps all data into one graph 
for the tested range of viscosities, frequencies, magnetic fields and layer thicknesses. 

PACS numbers: 47.20.-k, 47.20.Gv; 47.35. +i, 47.65. +a 



I. INTRODUCTION 

The Faraday instability denotes the parametric gen- 
eration of standing waves on the free surface of a fluid 
subjected to vertical vibrations. The study of this phe- 
nomenon dates back to the observations by Faraday in 
1831 Q. The initially flat free surface of the fluid be- 
comes unstable at a certain intensity of the vertical vibra- 
tions of the whole system. As a result of the instability, 
a pattern of standing waves is formed at the fluid sur- 
face. The typical response is subharmonic, i.e., the wave 
frequency is half the frequency of the excitation. The 
harmonic response can be observed on a shallow fluid 
at low frequencies Faraday waves allows one to in- 
vestigate symmetry breaking phenomena in a spatially 
extended nonlinear system. Therefore they experience a 
renewed interest in recent years. A detailed experimental 
study of the various patterns on a viscous fluid has been 
performed by Edwards and Fauve || who used a one- 
frequency as well as a two-frequency forcing. Parallel 
rolls, hexagons, and a twelvefold quasi-pattern were ob- 
served. Though typical Faraday waves are subharmonic, 
the twelvefold quasi-pattern appears to be harmonic in 
a certain region of parameters. Binks et al. have shown 
experimentally that the depth of the layer Q| and the 
excitation frequency JH) affect significantly the regions, 
where cither rolls or hexagons or squares are observed. 
In [|[ @, |[ |§] it is revealed that a two-frequency excita- 
tion leads to a great variety of wave patterns. In par- 
ticular, the observed patterns were triangles Jq] , super- 
lattices formed by small and large hexagons J/J, squares 
H 0, ||, ^| , and rhomboid pattern j^] . 

The comprehensive linear stability analysis of the Fara- 
day instability on an arbitrarily deep layer of a viscous 
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nonmagnetic fluid has been performed by Kumar and 
Tuckerman p0[ . This analysis was tested experimen- 
tally in flTT| and an excellent agreement between the pre- 
dicted and experimental data was found. Wcizhong and 
Rongjue |l2j extended the linear analysis pi for the case 
of an arbitrary periodic excitation. In j^, [Tj| the low fre- 
quency region is studied in great detail. Bicritical points, 
where transitions from one type of response to others oc- 
cur, are predicted and experimentally confirmed in [Q. 
In JlJ, |l5| an analogy between the Faraday instability 
and a periodically driven version of the Rayleigh- Taylor 
instability is exploited. In jl4| a scaling law is suggested, 
which satisfactorily describes the behavior of the system 
in a wide range of parameters. Kumar [ fl5| discusses the 
mechanism of the wave number selection in the Faraday 
instability on high-viscous fluids. 

In order to solve the pattern selection problem for 
the Faraday instability it is necessary to take into ac- 
count nonlinear interactions between different excited 
modes. The nonlinear behavior of the system has been 
studied theoretically in a great number of papers (see 
pq, [It], [j^, [l9[ and references therein). An ampli- 
tude equation for an infinitely deep layer of a viscous 
fluid is derived by Chen and Vinals |l9). A good agree- 
ment between the predicted frequencies for the transition 
between regions of different symmetry and the exper- 
imental observations |^] is found. 

Magnetic fluids (or ferrofluids) are colloidal dispersions 
of single domain nanoparticles in a carrier liquid. The at- 
tractiveness of ferrofluids stems from the combination of 
a normal liquid behavior with the sensitivity to magnetic 
fields. This enables the use of magnetic fields to control 
the flow of the fluid, giving rise to a great variety of new 
phenomena and to numerous technical applications |2l[| . 

One of the most interesting phenomena of pattern for- 
mation in ferrofluids is the Rosensweig instability j22) . At 
a certain intensity of the normal magnetic field the ini- 
tially flat surface of a horizontal ferrofluid layer becomes 
unstable. Peaks appear at the fluid surface, which typi- 
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cally form a static hexagonal pattern at the final stage of 
the pattern forming process p3| . By including vertical vi- 
brations to that setup, Miiller |Q analyzed the Faraday 
instability on viscous ferrofluids subjected to a vertical 
magnetic field in the frame of a linear stability analy- 
sis. It has been found that the joint action of the two 
destabilizing factors leads to a delay of the Rosensweig 
instability. By an appropriate choice of the system pa- 
rameters one can observe either a normal or anomalous 
dispersion. The predictions of p4| were confirmed exper- 
imentally in |25| . 

Parametric waves on the surface of a ferrofluid can be 
excited using a vertical |2(| or a horizontal |27], |28| al- 
ternating magnetic field. This phenomenon is called the 
magnetic Faraday instability. In p6| the dispersion rela- 
tion was investigated experimentally and it displays two 
significant features. The response of the surface waves is 
harmonic with respect to the frequency of the magnetic 
field and the wave vector of the resulting rolls is parallel 
to the field, i.e. the crests and troughs of the rolls are 
perpendicular to the field. In p8| a supercritical transi- 
tion from rolls to rectangles was observed and explained 
by means of a weakly- nonlinear analysis. 

In the present paper the stability of the surface of a fer- 
rofluid subjected to vertical vibrations and a static hor- 
izontal magnetic field is studied in a wide range of the 
system parameters. The horizontal magnetic field tends 
to decrease the curvature of the ferrofluid surface along 
the direction of the field (^9), i.e. it tends to stabilize the 
flat surface. On the other hand, vertical vibrations tends 
to destabilize a flat surface. Thus the aim of the linear 
stability analysis is to study the behavior of a ferrofluid 
subjected to two, in a certain sense competing factors. 
It will be shown that the magnetic field allows one to 
control the stability of the surface significantly and to af- 
fect the symmetry of the linearly most unstable pattern. 
This opportunity can be used in numerous technical ap- 
plications, where magnetic fields tangential to the fluid 
surface and perpendicular vibrations are typical. 



all three media 

divB«=0, rotH^^O, i= 1,2,3, (2.1) 

where superscripts denote the different media: 1 air, 
2 magnetic fluid, and 3 container. A linear relation 
between the magnetization of the fluid M' 2 ) and the 
strength of the magnetic field inside is assumed, M^ 2 ) = 

mo(/v - i)h( 2 ). 

The fluid motion is governed by the continuity equa- 
tion and the Navier-Stokes equations 

divv = 0, (2.2a) 
<9v 1 

— + (v grad) v = — gradp + i^Av + g(i) . (2.2b) 
ot p 

Here v = (u, v, w) is the fluid velocity, p is the pressure, 
and v is the kinematic viscosity of the fluid. The vertical 
vibrations add a periodic term to the gravity acceleration 
go, i.e., a modulated value g(t) = (0, 0, — g — a cos(wi)) 
appears in the Navier-Stokes equations. Here a is the 
acceleration amplitude and u is the angular frequency of 
the vibrations. The components of the stress tensor T^ 2 ' 
read 

T§ ] = I -P - Mo jMdH' - no^- I % (2 3) 

+ HiBj + pv(d l Vj + djVi) . 

The governing set of equations has to be supplemented 
by the boundary conditions. For the magnetic field the 
conditions are the decay of all perturbations far from the 
ferrofluid (z —> ±oo) and the continuity of the normal 
(tangential) component of the induction (strength) of the 
magnetic field across the air-fluid interface (z = Q and 
at the bottom of the container (z = —h) 



II. SYSTEM AND BASIC EQUATIONS 

A dielectric, viscous, and incompressible magnetic fluid 
with constant density p and permeability p, r is consid- 
ered. The laterally infinite ferrofluid layer of arbitrary 
depth h is subjected to a homogeneous dc horizontal mag- 
netic field and harmonic vertical vibrations (Fig. ^). The 
plane z — coincides with the nondeformed surface of 
the ferrofluid. The fluid layer is bounded from below by 
the bottom of the nonmagnetic container and has a free 
surface described by ((t, x, y) with air above. A homoge- 
neous magnetic field Ho = (Hq,0, 0) is applied along the 
x-axis. Due to zero electrical conductivity of the fluid, 
the static form of the Maxwell equations is used for the 
strength H and the induction B of the magnetic field in 




FIG. 1: A horizontally unbounded ferrofluid layer (2) is placed 
in a nonmagnetic container (3) with air (1) above. The system 
is subjected to a horizontal magnetic field Ho and harmonic 
vertical vibrations a(t). 
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H (1) = H at z^oo, (2.4a) 

B<p=Bg\ hW=hV> at z = (, (2.4b) 

B^=B^\ hW=hW at z = -h, (2.4c) 

H (3) = H at z^-oo. (2.4d) 

The subscripts n and r denote the normal and tangential 
components of the vector. 

The hydrodynamic equations are closed by the no-slip 
condition at the bottom of the container, 



at z 



(2.5) 



and the kinematic boundary condition at the free surface 
of the fluid 



w = dtC + (v grad)£ at z = ( 



(2.6) 



The equations for magnetic field and the fluid motion are 
coupled by the continuity condition for the stress tensor 
across the air-fluid interface, which completes the state- 
ment of the problem, 



{h ^ 
- p {1) + p + /io J MdH' + p ^- \ 
o ' 

—pvrii(diVj + djVi) — oKrij . 

(2.7) 

Here a is the surface tension, K — divn is the surface 
curvature, p^ is the atmospheric pressure above the fluid 
(assumed to be constant), and n is the unit vector normal 
to the surface given by 



grad[z - C(t,x,y)] 
|grad[z - ((t,x,y)}\ 



(2.8) 



III. LINEAR STABILITY ANALYSIS 

Following the standard procedure Jl^, |24|, the gov- 
erning equations and the boundary conditions have been 
linearized in the vicinity of the nonperturbed state 



v = 0, C = 0, H»=H , 
Po = P (1) - y M ° H ° ~ g ^ z ' 



1,2,3 , 



where M = (fi r — 1)Hq and po is the pressure in the un- 
perturbed state. The linearized governing equations for 
the small perturbations, which encompass the magnetic 



field strength Hi, the pressure pi, and a nonzero velocity 
v of the fluid, read 



jJLojj, r div H^ = 0, 



divv = 



1 



dv 

-rrr = — gradpi 
ot p 



rotH^ = 0, f = 1,2,3 , 

(3.1a) 
(3.1b) 

i/Av . (3.1c) 



Here the linear relation between the induction and the 
strength of the magnetic field is used. At the first order 
with respect to the perturbations, the boundary condi- 
tions are 







H« 


= at z — ► oo , 


(3.2a) 


H ln - 


f*r n ln ' 




= H$ at z = , 


(3.2b) 


fJ-r-Hin 


o-(3) 
- H ln > 




= h[^} at z = —h , 


(3.2c) 








= at z — > — oo , 


(3. 2d) 






v 


= at z = —h , 


(3.2e) 




w - dtC 


= at z = , 


(3.2f) 




d z u 


+ d x w 


= at z = , 


(3.2g) 




d z v 


+ 8yW 


= at z = , 


(3.2b.) 


g(t)C - 


-Pi- HqMqHx 


+ 2pvd z w = crAj_£ at z 


= 0, 










(3.2i) 


where A 


_L ®xx 


+ dyy. 


The stability of the flat surface 



with respect to standing waves is analyzed by using the 
Floquet ansatz for the surface deformations and the z- 
component of the velocity 



C(t,x,y) = sm(kr) e (*+ ia ")* £ Q n e ln ^ , (3.3a) 

n— — oo 
oo 

(t,x,y,z) = sin(kr)e( s+M ^ t Wn(z)e lnuJt , (3.3b) 



w 



n— — oo 



where k = (k x ,k y ) is the two-dimensional wave vector, 
s is the growth rate, and a is the parameter determin- 
ing the type of the response. For a — the response is 
harmonic whereas for a = 1/2 it is subharmonic. Expan- 



sions similar to (B.3b) are made for all other small per- 
turbations and are inserted into the linearized governing 
equations (3.1). The functions of the vertical coordinate 
in the Floquet expansion are given by linear combina- 
tions of the exponential functions e ±kz and e ±9 ™ z with 
= k 2 + [s + i(a + n)ui]/is and 8t(q n ) > 0. The condition 
of reality for £(i, x, y) leads to the equations pi 



C-ri Cn! 
C-n = Cn- 



a = , 
a = 1/2 



(3.4a) 
(3.4b) 



The boundary conditions (3.2) allow one to express all 



the perturbed quantities in terms of the coefficients 
which satisfy the equation 
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(W n ( n - aCn-i - e ^ a+n ^ t = , (3.5) 

n=— oo 

where 
W n = -2 



a = \F(a,k,u>,B eS , a, p, fi r )\ 



(3.10) 



k[q n coth(q n h) — k coth(kh)] 
(g„[4fc 4 + (k 2 + ql) 2 } coth(fc/i) coth(q n h) 

Aq n k 2 (k 2 + q 2 n ) 



2\2l 



k[Ak 2 q 2 n + (k 2 + q 2 n ) 



9o 



smh(kh) smh(q n h) 
ok 2 k(n r - l) 2 A(kh) fk x B ^ 1 



and 



A{kh) 



e kh (jj. r + l) + e- kh (iJ. r -l) 

(jW + l) 2 - e- kh (n r - 1) 



-J A;//. 



(3.6) 



(3.7) 



Here Bo = /^oHo is the induction of the applied magnetic 
field. 



The essential differences of Eqs. )— (3.7) in compari- 
son with the dispersion relation for the surface instability 
in a vertical magnetic field J24|, |3lJ are the following: 
First, the term proportional to B 2 has the same sign as 
the terms related to surface tension and gravity. This im- 
plies that a horizontal magnetic field alone cannot induce 
any instability. Second in the same term, the factor k x 
appears instead of k in the case of the normal magnetic 
field. This allows one to introduce the effective field 



B 



off : 



(kgp) 

k 2 



B cos(0)- 



(3.8) 



where 9 is the angle between k and B . Hence, the in- 
fluence of a magnetic field with induction Bo on a wave 
propagating along an arbitrary direction k is equivalent 
to the influence of a field with induction B c ff, which is 
parallel to k. 



Equation (3.5) has to be satisfied for all times which 
implies that each term of the sum equals to zero. Us- 
ing the re lati ons between £ n with positive and negative 
numbers (3^), one gets the set of equations 



WoCo - < 


- aCi = 0, 


a = , 


(3.9a) 


W ( - <o 


- aCi = 0, 


a =1/2, 


(3.9b) 


l — aCn-i — 


aCn+i = 0, 


n = 1, • • • oo . 


(3.9c) 



A cutoff at n = N (in the present work N — 100) leads to 
a self-consistent equation for the acceleration amplitude 



where F is a complex functi on ex pressed in terms of con- 
tinued fractions. Equation ( [3.10 ) can be solved numeri- 
cally and gives the dependence of a on A; at fixed param- 
eters. The critical values of the acceleration amplitude 
a c and the wave number k c correspond to an absolute 
minimum of the curve a(k) at zero growth rate (s = 0). 



IV. 



RESULTS AND DISCUSSION 



In the following the effective field is given in units of 
the critical induction B c r for the Rosens weig instability 
on an infinitely deep layer of ferrofluid |22j]. Figure |2j 
presents marginal stability curves for a viscous ferrofluid 
at low frequency. The dependence of acceleration am- 
plitude on the wave number for s — divides the phase 
space into regions, where the surface of the ferrofluid is 
stable or unstable with respect to parametrically driven 
standing waves. The principal data which can be ex- 
tracted are the critical acceleration amplitude, the wave 
number, and the number of the tongue to which they be- 
long. The number of a tongue I (from left to right) is the 
order of response: the basic wave frequency related to the 
l-th tongue is — lu/2. The odd and even tongues are 
the regions, where either a subharmonic or a harmonic 
instability develops. It can be seen that all tongues shift 
towards smaller wave numbers under the influence of an 
applied magnetic field. In the case presented in Fig. || 
the field also causes a transition from a subharmonic to 
a harmonic response. 

The dependencies of the critical acceleration amplitude 
and the critical wave number on the excitation frequency 
f = lu/2tt are presented in Fig. ^|. In the high frequency 
region, the instability is subharmonic, i.e. 1 = 1. In the 
small frequency region bicritical points appear, which is 
why the dependence of a c (f) is not smooth and the de- 
pendence of k c (f) is discontinuous (see inserts in Fig. |^). 
Bicritical points are those points in the parameter space, 
where the absolute minima of a(k) is equal to the local 
minima of two neighbouring tongues. In the zero field 
case for instance, such an overlap happens between the 
first subharmonic tongue (1 = 1) and the first harmonic 
tongue (I = 2) at /b c ,i — 1-86 Hz. For frequencies below 
/bc,i, the tongue of 2nd order gives the critical acceler- 
ation until the next bicritical point at /b c ,2 — 1-55 Hz. 
Below /be. 2 the tongue of 3rd order gives there the lowest 
threshold until the third bicritical point and so on. Thus 
with decreasing frequency subsequent transitions from a 
tongue of the order I to a tongue of the order I + 1 occur. 

It is seen in the insert in Fig. ||(b) that the magnetic 
field shifts the bicritical points. If the fluid viscosity is 
low (as in Fig. |^), then the magnetic field decreases the 
frequencies of the transitions. For higher viscosities, e.g., 
v = 10 -4 m 2 /s, the bicritical points are shifted towards 



higher frequencies. This is exemplarily shown in Fig. ||for 
the frequency of the first bicritical point, /b c ,i> for a high- 
viscous fluid (curve 1) and a low-viscous fluid (curve 2). 
For frequencies above (below) the curves the response of 
the system is subharmonic (harmonic). Thus for a high- 
viscous fluid and an excitation frequency of / = 10 Hz 
at zero applied field a subharmonic response (CI = u>/2) 
is expected, whereas at B c ff = B c r a harmonic response 
(CI = lj) precedes. For ordinary (nonmagnetic) fluids 
bicritical points were also found in the case of both low- 
viscous fluids (2| and high- viscous fluids p3[ . 

A feature of the dependence of a c (f) [Fig. ||(a)] is the 
appearance of a pronounced minimum for all tested fields. 
To explain this behavior, one has to recall that viscous 
damping is the reason for the finite (nonzero) value of 
the critical acceleration amplitude. The damping occurs 
due to the stresses in the bottom layer and in the bulk 
fluid. The dimensionless quantity d = k c h can be used 
to determine which part of the damping is predominant. 
For a shallow fluid layer the inequality d < 1 holds and 
therefore the damping in the bottom layer is predomi- 
nant (first regime). For a deep layer, where the relation 
d 3> 1 is fulfilled, the damping in the bottom layer is 
of no importance and the damping in the bulk fluid is 
predominant (second regime). 

The first regime occurs at the low frequency region 
in Fig. ||, where k c < 500 m _1 . Inside this region, as 
the frequency increases the effective depth d of the fluid 
is increasing, too. Therefore the viscous stress in the 
bottom layer is weakening and consequently the critical 
acceleration amplitude decreases. The second damping 
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FIG. 2: Neutral stability curves for the excitation frequency 
/ = 10 Hz and layer depth h — 2 mm. Dashed (solid) lines 
correspond to B c s = (B cR ). Filled and unfilled tongues 
represent regions of subharmonic and harmonic responses. a c 
and k c (a' c and k' c ) are the critical acceleration amplitude 
and critical wave number for B a g — B cR (B c a = 0). The 
parameters of the fluid are v = 10 -4 m 2 /s, a = 0.0265 N/m, 
p = 1020 kg/m 3 , n r = 1.85, and B cR = 17.28 mT. 



30 
_ 20 
^10 




N J 20 


B 










1, 

1 1 

\\ 1 












D 








1.5 


2 
















\vc 












D 










(a) 



20 40 60 80 100 

/(Hz) 



1500 



1000 



500 




40 60 

/(Hz) 



100 



FIG. 3: Frequency dependencies of the critical acceleration 
amplitude a c (a) and the critical wave number k c (b) for h = 2 
mm. Inserts: low-frequency behavior of the quantities. The 
effective field is B = 1.5B cR (curve A), B cR (B), 0.5B cR (G), 
(D). The fluid parameters are fx r = 1.85, a = 0.0265 N/m, 
v = 5.88 x 10~ 6 m 2 /s, p = 1020 kg/m 3 , and B cR = 17.28 mT. 



regime is typical for waves with k c above 500 m -1 . Since 
dissipation in the bulk fluid is proportional to k\ [ p2| , 
the damping becomes stronger and as a result the criti- 
cal acceleration amplitude increases with frequency. The 
transition from the first to the second damping regime 
leads to the nonmonotonic dependence of the critical ac- 
celeration amplitude on the excitation frequency. 

The most pronounced effect caused by the magnetic 
field is the decrease of the critical wave number [cf. curves 
D-A in Fig. ||(b)] as already observed in [^9|. The re- 
duction of the critical wave number affects the critical 
acceleration amplitude. This influence is different at low 
and high frequencies due to the following reason. The 
decrease of k c results in (i) the reduction of the effective 
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FIG. 4: Frequency of the first bicritical point /bc,i versus the 
effective magnetic induction. For a high-viscous fluid (y = 
10 -4 m 2 /s, curve 1) an increasing induction increases /bc,i, 
whereas for a low- viscous fluid [y = 5.88 x 10" 6 m 2 /s, curve 
2) /bc,i is decreasing. 



fluid depth d, which intensifies the stress in the bottom 
layer and (ii) the weakening of the dissipation in the bulk 
fluid. Therefore in the case of the first damping regime 
the viscous dissipation is intensified, whereas in the case 
of the second regime it is reduced by the magnetic field. 

In the low frequency region (d = k c h < 1), the first 
damping regime occurs independent of the magnetic field. 
Therefore the critical acceleration amplitude increases 
with the growth of the effective field. 

At high frequencies the dependence of the critical accel- 
eration amplitude on the effective field is more intricate. 
In Fig. |5|(a) the dependence of a c (B e g) is presented for 
a frequency / = 100 Hz. It is seen that for an infinitely 
deep fluid the critical acceleration amplitude decreases 
monotonicly with the increase of the magnetic field. In 
the case of a finite depth of the layer, there are two dif- 
ferent forms of the dependence of a c on B e g. In the case 
of h > 1.14 mm a minimum in the critical acceleration 
amplitude is observed. It implies that a moderate mag- 
netic field lowers the threshold value of the acceleration 
amplitude whereas a strong magnetic field stabilizes the 
surface, i.e., the critical acceleration amplitude is increas- 
ing. With the decrease of the layer depth [from curve D 
to A in Fig. ^(a)], the observed minimum shifts to the 
lower fields and becomes less pronounced. At h < 1.13 
mm the critical acceleration amplitude monotonicly in- 
creases with the effective field. 

To explain these qualitative changes in the dependence 
of a c (B c s) it is useful to consider d(B c g) [see Fig. |(b)]. 
The analysis shows that in the region d ^> 1 (the second 
damping regime) the decrease of the critical wave num- 
ber caused by the magnetic field leads to a weakening 
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FIG. 5: Critical acceleration amplitude a c (a) and the dimen- 
sionless depth d — k c h of the fluid (b) versus the induction of 
magnetic field for / = 100 Hz. The layer depth is 1 mm (curve 
A), 1.5 mm (B), 2 mm (CT), 5 mm (D), and the curve (E) is 
obtained for an infinitely deep layer. The fluid parameters 
are those of Fig. ^. 



of the viscous damping. Therefore the critical acceler- 
ation amplitude is decreasing with the increase of the 
field until the effective depth d of the fluid becomes of 
the order of unity. A further decrease of d caused by 
the magnetic field increases the viscous damping in the 
bottom layer and consequently the critical acceleration 
amplitude. Thus, the transition from the second damp- 
ing regime to the first one results in the nonmonotonic 
dependence of a c (B g) [Fig. ||(a)]. The effective field, 
where the transition occurs, is decreasing with decrease 
ofh, and in the case of a very shallow fluid (here, h < 1.13 
mm), the second damping regime is not observed at all. 

The nonmonotonic dependence of a c (B c ff) allows one 
to select the type of the linearly most unstable pattern 
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(corresponding to a minimal a c ) by applying a horizon- 
tal magnetic field. When the induction B of the ap- 
plied field is smaller than B* s [a c (B e g] is minimal at 
B* s ] Faraday waves with k| |Bp are favorable. In this 
case B e ff = Bo [see Eq. (3.8)], i.e., B e s has a maximal 
possible value at given Bq- It implies that the first un- 
stable pattern are rolls perpendicular to the magnetic 
field. If Bq > B* s then the angle between the favorable 
wave vector of perturbation and the applied magnetic 
field satisfies the relation cos6> = ±B* S /B . Thus, for 
B > B* s ^ rolls with an angle +8 or an angle —9 or a 
rhombic pattern as a superposition of both is most prob- 
able. If the layer depth is small enough B* s becomes zero 
[see curve A in Fig. |(a)]. That is equivalent to cos6> = 0, 
i.e., rolls parallel to the magnetic field are expected. 

The parametrical generation of surface waves by a hori- 
zontal alternating magnetic field was observed in |^7|, . 
In these studies the driving force is the magnetic field and 
it is anisotropic. Only perturbations with a wave vector 
parallel to the field can be unstable in the linear approx- 
imation. Therefore the first unstable pattern observed in 
the experiments were always rolls perpendicular to the 
field. In Q a supercritical transition was observed from 
rolls to a rectangular pattern. This transition is caused 
by nonlinear interactions of the different modes. In the 
present study there is no anisotropic driving and there- 
fore Faraday waves along an arbitrary direction can be 
excited. As a consequence, different patterns can be gen- 
erated (rolls along an arbitrary direction or a rhombic 
pattern as discussed above). 

In the case of low viscosity the obtained results for 
the critical acceleration amplitude are in the very good 
agreement with an approximation suggested by Muller 
et al. (Eq. (9) in §] and Eq. (4.1) in @). It should 
be noted that the influence of a normal magnetic field 
(studied in |24|, |26|) and a horizontal field (studied in the 
present paper) on Faraday waves are different. The nor- 
mal magnetic field increases the critical wave number of 
the Faraday waves ]2(| , whereas the horizontal magnetic 
field decreases the wave number. Thus, the sensitivity of 
fcrrofluids to magnetic fields allows one both to increase 
and to decrease the critical wave number. In this way the 
relative importance of the stress in the bottom layer and 
the dissipation in the bulk fluid can be changed. There- 
fore magnetic fields are a convenient way to control the 
stability of the surface. 

The Faraday instability on nonmagnetic fluids with 
high viscosity has been investigated in Q . The authors 
suggest a scale for the acceleration based on an analogy 
between the Faraday instability and a periodically driven 
version of the Rayleigh- Taylor instability. They observe 
a data collapse in the range of the dissipation parameter 
(S/h) 2 from 0.1 to 0.3, where S = \fvjbj is the dissipative 
length scale. In this parameter range our results are in 
a good agreement with the suggested scaling law. How- 
ever, at low frequencies the scaling behavior suggested 
in [jl4| strongly underestimates a c . Additionally, the re- 
maining parameters neglected in jl4| become essential at 
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FIG. 6: Scaled critical acceleration amplitude versus the di- 
mensionless depth of the layer. Data are plotted for seven 
different sets with various combinations of the fluid viscosity, 
excitation frequency, effective field, and the layer thickness 
(see Table . The remaining parameters are those of Fig. |§| 



TABLE I: Sets of the parameters for data plotted in Fig. 



Set 


V 


(m 2 /s) 


/ (Hz) 


Bcff {B c r) 


h (mm) 


Symbol 


1 


5.8 


8 x 10 -6 


5 


1 


0.5 .. . 600 


+ 


2 




10 -3 


10 





5 . . . 100 


V 


3 




10" 4 


100 





1 . . . 100 


□ 


4 


5.8 


8 x 10 -6 


100 


... 3 


2 


A 


5 


5.8 


8 x 10~ 6 


1 . . . 100 


1 


5 


< 


6 




10~ 3 


5 . . . 100 


1 


2 




7 


10 _ 


6 ...HT 3 


48 


1 


2 


o 



kch < 1. Therefore, the scaling propagated in [|14j is of 
limited use for ferrofluids in magnetic fields. 

The aim now is to reveal a scaling law for the thresh- 
old values of the Faraday instability. The behavior of 
the critical acceleration amplitude as a function of the 
parameters of the system is determined by the dimen- 
sionless depth d of the fluid (see Figs. ||, ||). Therefore, it 
is natural to look for a scaling law of the form a c = d c (d), 
where d c is the scaled acceleration. As h — > oo the re- 
quired function should approach an asymptotic value. 
Thus, the appropriate scale for the acceleration ampli- 
tude is the value a c ,oo of the critical acceleration for an 
infinitely deep fluid. The quantity k c flv is a reasonable 
estimation for a c ,oo |l4) . Introducing this scale for the 
acceleration amplitude d c = a c / '(vQk c ), it is possible to 
map all the above presented dependencies in a single plot 
as it has been done in Fig. 0. 

Figure |^ presents the dimensionless amplitude as a 
function of the effective depth d = k c h of the fluid. The 
latter can be varied by means of changing the physical 
depth of the layer (pluses, down triangles, and squares 
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in Fig. ^]), varying the applied magnetic field (up trian- 
gles), excitation frequency (left and right triangles), and 
the viscosity of the fluid (diamonds). It is seen that for 
a wide range of parameters the dependencies of a c re- 
lated to different fluids and varying parameters are in a 
rather good agreement with each other. The only devi- 
ation from the common behavior appears for a low fre- 
quency (pluses) at intermediate d. For d up to unity the 
dimensionless acceleration amplitude can be reasonable 
approximated by d c w 2ld~ 2 6 (solid line in Fig. ^). For 
d > 1 an estimation a coo ks 7 approximates the exact 
results with a maximal error of about 12 %. It is worth 
noting that the dependencies of the critical acceleration 
amplitude on the surface tension and the fluid density 
follow the same common behavior as shown in Fig. ||. 

V. CONCLUSION 

The linear analysis of the Faraday instability on a vis- 
cous ferrofluid subjected to a horizontal magnetic field 
has been performed. A horizontally unbounded ferrofluid 
layer of a finite depth has been considered. The depen- 



dencies of the critical acceleration amplitude a c and the 
critical wave vector on the excitation frequency / and the 
induction B c g of the magnetic field have been obtained 
for different depths of the layer in a wide range of fluid 
viscosities. The regions have been found, where the vis- 
cous stress either in the bottom layer (the first regime) 
or in the bulk fluid (second regime) are predominant. A 
transition from the second damping regime to the first 
one can be caused by decreasing the excitation frequency 
or by applying a horizontal magnetic field. The transi- 
tion results in the nonmonotonic dependencies of a c (f) 
and a c (B c g). It is shown that one can select the first 
unstable pattern of Faraday waves by the appropriate 
choice of the depth of the ferrofluid and the induction 
of the magnetic field. A scaling law is suggested, which 
describes the behavior of the system in a wide range of 
the system parameters. 
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